Driving-induced bistability in coupled chaotic attractors 
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We examine the effects of symmetry-preserving and breaking interactions in a drive-response 
system where the response has an invariant symmetry in the absence of the drive. Subsequent to 
the onset of generahzed synchronization, we find that there can be more than one stable attractor. 
Numerical, as well as analytical results establish the presence of phase synchrony in such coexisting 
attractors. These results are robust to external noise. 
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I. INTRODUCTION 

Synchronization is an ubiquitous phenomenon in the 
dynamics of interacting nonlinear systems. Although 
first used to denote the establishment of identical dy- 
namics in weakly coupled systems, in recent decades it 
has been realized that the concept extends itself to de- 
scribe correlated motion in coupled systems quite natu- 
rally. Based on the nature of the systems and the manner 
in which they are coupled, several types of synchroniza- 
tion are known [i|-|8|. 

Systems can be coupled either mutually or one-way as 
in the "drive-response" scenario wherein the response is 
enslaved by the forcing (master) system. In this latter 
case, the equations of motion can usually be represented 
in the following skew-product form. 



X = F(X) 



G(Y,X), 



(1) 



(2) 



where X and Y are the dynamical variables of the drive 
and the response systems respectively, and F and G are 
continuous and differentiable vector fields. The occur- 
rence of generalized synchronization (GS) in such sys- 
tems implies the existence of a unique (but not neces- 
sarily differentiable) functional relationship Y = ^(X) 
between the variables of the drive and response systems 
[Sl . The onset of GS occurs when the largest conditional 
Lyapunov exponent (LCLE) of response system becomes 
negative [9|,|l0|- The presence of GS is further determined 
by constructing an auxiliary copy (Y^) of the response 
system as 



Q . This form of synchrony is quite common in natural 
systems wherein there is a range of time-scales in the 
dynamics. 

Studies of GS have largely focused on chaotic sys- 
tems having single-scroll dynamics where the oscillation 
is around a single fixed point. However systems with 
multi-scroll dynamics when the trajectory hops around 
more than one fixed point are also known, and are also 
known to be shown by the important class of systems 
having inherent symmetries. Guan et a/. [12] have consid- 
ered such symmetric systems in drive-response coupling 
scenario and observed bistability within the GS regime. 
However, the dynamics of these bistable attractors has 
not been studied in detail. 

The motivation of the present work is the study of in- 
duced bistability in unidirectionally coupled symmetric 
response systems. The coupling may either preserve or 
break the symmetry of the response system, but we find 
that the resulting bistable attractors are always phase- 
synchronized. Our results further elucidate the nature of 
GS and extend the auxiliary system approach to the anal- 
ysis of phase synchrony between response system and its 
auxiliary copy. These observations are seen to be robust 
to noise. 

The paper is organized in following manner. The possi- 
ble coupling schemes in context of symmetry of response 
system are discussed in Sec. II. This is followed by nu- 
merical and analytical demonstration of phase relations 
with different coupling schemes in Sec. HI and IV. The 
robustness of these results is verified in Sec. V. The paper 
concludes with a summary in Sec. VI. 



Y' = G(Y',X), 



(3) 



where the occurrence of GS between drive (X) and re- 
sponse (Y) is confirmed by observing complete synchro- 
nization (CS) between the response system (Y) and its 
auxiliary copy (Y') |11| . 

When non-identical systems are coupled, they can ex- 
hibit a form of correlated dynamics that is termed phase 
synchronization (PS): the phases of oscillation can be- 
come entrained while the amplitudes remain uncorrelated 



II. COUPLING SCHEMES 



In present work, we are interested in studying the re- 
sponse behavior of a particular class of systems, namely 
those which are invariant with respect to some symmetry 
transformations. A familiar example is provided by the 
Lorenz oscillator [13| , the governing equations of which 



are 



X = a{y — x) 
y = rx — y — xz 

z = xy — j3z. 



(4) 



This system is invariant under the transformation T : 
{—x^—y^z) -^ {x^y^z), and thus the attractor is sym- 
metric about the z- axis, or equivalent ly it has inversion 
symmetry in the x — y plane. As is well known, this sys- 
tem has a "butterfly" shaped strange attractor (at pa- 
rameters cr=10, r=28, and /3=8/3), where the trajectory 
moves between quadrants 1 and 3 in x — y plane. 

We consider driving this system by an external dy- 
namical system, with diffusive scalar coupling [1^ , where 
the drive is coupled to a single variable of response sys- 
tem. Clearly the response system can be forced in a 
manner that either preserves or breaks the symmetry. In 
the present case, if the (unidirectional) coupling is intro- 
duced in the z variable then the symmetry in x and y 
is preserved, and this coupling is termed symmetry pre- 
serving. If the response is coupled through either the x 
or the y variables, the inversion symmetry with respect 
to the z axis or x — y plane is lost and this is a symmetry 
breaking coupling. The dynamical behavior of response 
system under these two coupling schemes are studied in 
next sections. 



III. SYMMETRY-PRESERVING 
INTERACTION 

Consider the Lorenz system (Eq. |4]) driven by a chaotic 
Rossler j]L5| oscillator, the complete system being given 
by 



xi = -yi - zi 

yi = xi + ayi 

ii = 6 + zi{xi - c) 

X2 = cr{y2 -X2) 

1/2 = rx2 -y2- X2Z2 

Z2 = X2y2- Pz2^e{zi- Z2). 



(5) 



The symmetry of the response Lorenz system in the 
X2 — y2 plane is preserved by the coupling. Fixing the 
internal parameters of drive and response oscillators to 
standard values, a=0.2, 6=0.2, and c=5.7 and cr=10, 
r=28, and /3=8/3 respectively, both systems have chaotic 
dynamics. The largest conditional Lyapunov exponent 
(LCLE) [9|, [lO| as a function of the coupling parameter e 
is shown in Fig. [IJa). With the increase in coupling pa- 
rameter the LCLE becomes negative at e ?^ 0.767, which 
indicates the onset of GS between the drive and response 
systems. By constructing an auxiliary system [11| and 
measuring the average synchronization error (A) between 
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FIG. 1: Lorenz (response) system: (a) the largest conditional 
Lyapunov exponent (LCLE) and (b) the average synchro- 
nization error (A) between response system and its auxiliary 
copy, with coupling parameter e. 



the response system and its auxiliary copy confirms GS. 
The error A [16| is defined as 



A = {V{{X2 - X3)2 + (2/2 - 2/3)2 + (^2 _ ^3)2))^ 

where (•) denotes a time average as well as over an en- 
semble of initial conditions, typically taken to be 100. 
The value of A as a function of e are shown in Fig. mb), 
and at the transition to GS, A takes two values, indicat- 
ing that there are two distinct attr actors, one on which 
A = and one on which A > 0. The trajectories cor- 
responding to the different attr actors are shown in Figs. 
[2fa) and (b), for e=0.74 and 1.59 respectively. 

Figure[2fa) shows that there is always a single attractor 
irrespective of initial conditions, prior to the onset of GS, 
but as can be seen in Fig. [2fb), the response system con- 
sists of two symmetric coexisting attractors A and B that 
are spatially separated. Symmetry preservation is clearly 
shown in Figs. [2jc) and (d) corresponding to attractors 
A and B respectively. With different initial conditions, 
if both the response system and its copy evolve to either 
the attractor A or attractor B then A = 0. However it is 
also possible that the trajectories of the response and the 
auxiliary converge to different attractors, and this gives 
a positive average synchronization error, namely A > 0, 
even though there is GS. 

This is termed as symmetry induced bistability, and 
the occurrence of two coexisting attractors A and B can 
be verified for a number of different initial conditions. 
One question of interest is the phase relationship between 
the distinct attractors (when A > 0). Shown in Figs. 
[3l^a) and (b) are the trajectories in X2 — xs plane (to 
infer the relative phase) when both the response system 
and its copy, lead to the same attractor (here A) or to 
different attractors {A and B) respectively for two sets 
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FIG. 2: (a) Lorenz attractor before the onset of GS at cou- 
pling strength e=0.74, (b) The symmetric pair of coexisting 
attractors A and B in response system for model Eqs.([5j), (c) 
, and (d) the expanded views of coexisting attractors A and 
B. Coupling parameter e=1.59. 



of different initial conditions. Fig. [3l^a) shows that the 
response system and its auxiliary copy are in-phase or 
in complete synchrony (x2 = X3,?/2 = ^3) which gives 
A = 0, while, Fig. [3l^b) presents the case when different 
initial conditions settle onto coexisting attractors (A and 
B\ The response and the auxiliary are in anti-phase 
synchrony (APS) since X2 = —^3 and y2 = —^3 . 



FIG. 3: Lorenz attractor: the evolution of response system (x 
variable) for two different sets of initial conditions at e=1.59, 
where these conditions converges to (a) the same attractor B 
(CS, X2 = X3), (b) either of symmetric coexisting attractors 
A and B {APS, X2 = —xs). The parameters are that of the 
Fig. 2(b). 



A phase variable can be defined for the the Lorenz 
system as (j){t) = arctan(z — zo)/{u — uq) by considering 
the projection of trajectory onto the {u^ z) plane where 
u = Y^x^ + y'^ and {uq^ zq) is the enclosed fixed point \Vq . 
We denote the phases of the response and auxiliary as ^2 
and 03 respectively. To establish the phase relationships 
between the response and its copy we begin with the 
trigonometric identity, 



tan( 



5) 



tan( 



tan( 



1 + tan (j)2 tan ^3 
using which one can write the phase relation as 



(6) 



Z2 - Z02 Z3 - Zos 



tan(02 - h) 



U2 - U02 Us - Uqs 



1 



Z2 - Z02 Zs - Zos 



U2 - U02 Us - Uqs 
{us - Uos){z2 - ^02) - (^2 - U02){ZS - ^03) 



(7) 



{U2 - U02){US - Uos) + {Z2 - ^02)(^3 - ^03) ' 

For any A > 0, when (/)2 — ^3 is or tt, tan(02 — 03)=O, 
giving 

(^3 - Uos){z2 - Z02) - {U2 - U02){ZS - Zqs) = 0. (8) 

Since the Lorenz attractor has the inversion symmetry 
T : (—X, —7/, z) -^ (x, y^z) in x — y plane, the z variables 
of the response system and the auxiliary remain identical 
in GS regime. Setting zs = ^2 and 2:03 = 2:02 Eq. ([8]) can 
be rewritten as 



{Z2 - ^02) [(^3 - ^03) - (^2 - ^02)] = 0- 



(9) 



Here, Z2 7^ 2:02 because 2:02 is the z— component of the 
unstable fixed point. It is easy to show that this leads to 



12 




_(b) 


1 ' 1 ' ' 


1 1 


_ 


. > 


- 




v«5. 






^ 


1 , 1 


1 



FIG. 4: Lorenz (response) system in symmetry breaking case: 

(a) the largest conditional Lyapunov exponent (LCLE) and 

(b) average synchronization error (A), with coupling param- 
eter e. 



the condition 

{xs + X2)(x3 - X2) + (^3 + y2){y3 - ^2) = 0, (10) 

which can be satisfied for both CS, namely (/)2 — </>3 = 
when xs=X2 and ^3=^2 as well as for xs = —X2 and 
Us = —y2 namely APS. Both cases are realized by our 
numerical results. 

IV. SYMMETRY-BREAKING INTERACTION 
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FIG. 5: (a) Asymmetric pair of coexisting attractors C and 
D in (x- coupled) symmetry breaking case, (b) and (c) the 
expanded views of coexisting attractors C and D. Coupling 
parameter e = 2.0. 



Consider the coupling in x in the response, namely the 
equations of motion 



X2 = cr{y2 - X2) + e{xi 
7/2 = rx2 -y2- X2Z2 

Z2 = X2y2 - PZ2. 



X2) 



(11) 



The internal parameters of drive and response systems 
are kept the same as in Eq. (|5]). Figs.lU^a) andlU^b) show 
the LCLE and A respectively for this coupling scheme. 
As discussed in earlier section, the fluctuation of aver- 
age synchronization error (A) having the values A > 
implies the presence of bistable chaotic attractors. As 
shown in Fig. [5ja), the response system consists of two 
coexisting asymmetric attractors C and D (in GS regime) 
for two sets of different initial conditions. However, these 
attractors, C and D^ are asymmetric in phase space as 
clearly shown in expanded Figs. [5fb) and (c) respectively. 

The synchrony properties of the asymmetric coexisting 
attractors C and D when A > in the GS regime are 
of interest. Fig. [6l^a) is for the case when both response 
and auxiliary lead to the same attractor (denoted D) for 



two sets of different initial conditions: the response and 
auxiliary are in-phase or in complete synchrony {x2 = 
^3,^2 = ^3)- Fig. [6l^b) is for the case when the response 
and auxiliary converge to distinct attractors: the two 
systems are in-phase. 



V. PERSISTENCE WITH NOISE 

In order to see the constancy of results obtained in 
case of symmetry preserving and breaking interactions, 
the effects of added noise are considered. Thus in the 
symmetry preserving case (cf. Eq. ([5])) we consider the 
system 

X2 = cr{y2 -X2) +<fS'i 

y2 = rx2 -y2- X2Z2 + ^S2 

Z2 = X2y2- /^Z2^e{zi- Z2)^^S3, (12) 



(a) 




FIG. 6: Symmetry breaking interaction: the evolution of re- 
sponse system (x variable) for two different initial conditions, 
where these conditions converge to (a) the same attractor C 
(A = 0), (b) either of asymmetric coexisting attractors C and 
D (A > 0)). Coupling parameter e = 2.0. 



and for symmetry breaking case (cf. Eq. (pT]) ). 

X2 = cr{y2 - X2) ^ e{xi - X2) ^ ^S^ 

y2 = rx2 -y2- X2Z2 + ^Sb 

Z2 = X2y2 - PZ2 + (,Sq. (13) 

Here Si {i = 1,2,.., 6) is (5-correlated random noise of 
strength ^. Other parameters (a, r, /3) are of the same 
values as used earlier in Sections II and III. The effects 
of noise for symmetry preserving and symmetry break- 
ing interactions are illustrated in left and right panels in 
Fig. [71 respectively with noise strength, ^ = 10.0. The 
variation of LCLE with coupling parameter, e are shown 
in Figs. [71(a) and (b) which clearly show that there is no 
significant difference as compared to Figs. [11(a) and[ll(a) 
respectively. The fluctuations of average synchronization 
error, A (Figs. [3 (c) and (d)) infer the occurrence of 
driving induced bistability even in presence of noise. 

The phase relations between coexisting attractors are 
shown in Figs. [3 (e)-(h). Figs. [71(e) and (f) show the 
PS state in presence of noise, similar to Figs. [3l (a) and 
[6l (a) respectively. Figure [71 (g) indicates the APS state 
similar to Fig. [31(b), while Fig|71 (h) shows the PS state 
as in Fig. [6l (b). We have also obtained the similar re- 
sults even with higher noise strength. Thus the results 
obtained in section II and III are robust to added noise, 
and hence implying the possibility of theirs experimental 
verification. 



VI. SUMMARY 

In present work we have studied the effect of external 
forcing on systems with invariant symmetry, using the 
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FIG. 7: The left and right panels correspond to the sym- 
metry preserving and symmetry breaking interaction respec- 
tively with noise (^ = 10.0). (a) and (b) show the largest 
conditional Lyapunov exponents (LCLE) while (c) and (d) 
show the average synchronization error (A) as a function of 
coupling strength, e. (e) and (f): the trajectories in (x2, xs) 
plane, when the response and its auxiliary copy converge to 
the same attractors B and D respectively, (g) and (h): the 
trajectories in (x2, X3) plane, when the response and its aux- 
iliary copy converge to the different attractors A or B and C 
or D respectively. The coupling parameter in (e-h) is e = 2. 



Lorenz equations as a prototypical example. The cou- 
pling has been introduced in two different ways, to keep 
the inherent symmetry of the driven system, or to break 
the symmetry. 

The onset of generalized synchronization between the 
drive and the response is studied, and we find that there 
is bistability after GS has set in: this is clearly drive in- 
duced. The coexisting attractors have been found to be 
symmetric or asymmetric depending upon the coupling, 
and we use the auxiliary system approach to verify our 
results. In the symmetry preserving case, we find that the 
attractors in the response and the auxiliary can be either 
in phase or antiphase, and we are able to present analyti- 
cal arguments for these cases. However for the symmetry 



breaking interaction the coexisting attractors (which are 
asymmetric) can only have in-phase synchrony. 

In this paper we have only used the chaotic Rossler at- 
tractor as a drive, but it should be noted that results are 
similar for any drive that has suitable internal timescales. 
Further, although, we have taken the response to be the 
Lorenz system here, we have observed similar results for 
phase synchronization states in other inversion symmet- 
ric response systems [18], and we believe that our results 
will hold quite generally. 
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